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ABSTRACT 

The IMACS Cluster Building Survey uses the wide field spectroscopic capabilities of the IMACS 
spectrograph on the 6.5m Baade Telescope to survey the large-scale environment surrounding rich 
intermediate-redshift clusters of galaxies. The goal is to understand the processes which may be 
transforming star-forming field galaxies into quiescent cluster members as groups and individual 
galaxies fall into the cluster from the surrounding supercluster. This first paper describes the survey: 
the data taking and reduction methods. We provide new calibrations of star formation rates derived 
from optical and infrared spectroscopy and photometry. We demonstrate that there is a tight relation 
between the observed star formation rate per unit B luminosity, and the ratio of the extinctions of 
the stellar continuum and the optical emission lines. With this, we can obtain accurate extinction- 
corrected colors of galaxies. Using these colors as well as other spectral measures, we determine new 
criteria for the existence of ongoing and recent starbursts in galaxies. 
Subject headings: galaxies: clusters, evolution, star formation, cosmic evolution 



1. INTRODUCTION 

That the properties of galaxies differ with environ- 
ment has been recognized at least since Hubble (1936). 
That much of the difference between clusters and the 
field is of recent origin has been known since Butcher & 
Oemler (1978). However the mechanisms that have pro- 
duced these population differences are still in dispute. 
Many processes have been suggested which are capable 
of transforming field-like populations into the predomi- 
nantly early Hubble types seen in clusters today. These 
include gas stripping by galaxy-galaxy collisions (Spitzer 
& Baade 1951), gas stripping by ram pressure from the 
intracluster medium (Gunn & Gott 1972), a shutoff in 
gas replenishment (Larson, Tinsley, & Caldwell 1980), 
tidal shocks, either due to the cluster core (Byrd & Val- 
tonen 1990, Henriksen & Byrd 1996), to unvirialized sub- 
clusters (Gnedin 2003), or to other galaxies (Richstone 
& Malmuth 1983, Icke 1985, Moore et al. 1996), and 
galaxy-galaxy mergers (Dressier et al. 1999, van Dokkum 
et al. 1999). Since all these processes produce, by design, 
the same outcome, and have, again by design, a qualita- 
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tively similar dependence on environment, distinguishing 
between them is challenging at best. Given that most ex- 
tant observations consist of snapshots at one epoch of the 
cores of individual clusters, it is hardly surprising that 
the responsible process(es) have still not been unambigu- 
ously identified. 

The IMACS Cluster Building Survey is an attempt to 
resolve this problem by following the evolution of galaxies 
as they move from the supercluster environment, through 
the infall stage and end finally with incorporation into 
the virialized central cluster. It takes advantage of the 
very wide field and high multiplexing of the IMACS spec- 
trograph on the Baade Telescope (Dressier et al. 2011), 
which allow one to observe, in one exposure, hundreds of 
galaxies over a 30' field, equivalent to a circle with radius 
of about 5 Mpc surrounding a cluster at redshift 0.4. The 
goal is to identify the changes that occur in galaxies as 
they move from field-like environments into increasingly 
dense and massive assemblies. The much longer baselines 
of time and environment which these data provide should 
help distinguish between the various candidate processes 
for transforming galaxies. 

In this paper we describe the design and execution of 
the survey, up through the production of redshifts, lu- 
minosities, colors, masses, and star formation rates of 
galaxies. We describe new calibrations of star forma- 
tion rates based on optical and infrared photometry and 
spectroscopy and we discuss several methods for detect- 
ing starbursts using the available data. The immediately 
following papers (Dressier et al. 2013- Paper II, Oemler 
et al. 2013- Paper III, Gladders et al. 2013- Paper IV) 
discuss certain aspects of the behavior of galaxies in 
the cluster and supercluster environment, and the evo- 
lution of the field galaxy population. Two papers using 
ICBS data to analyze environmental variations in the 
mass function of galaxies have already been completed 
(Vulcani et al. 2012, 2013). Future papers will address 
other aspects of both the field and cluster populations. 
Throughout this and following papers we shall assume 
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cosmological parameters of H Q = 71, ^matter — 0.27, 
tttot = 1.0. 

2. OBSERVATIONS 

We wish to map a group of rich intermediate rcdshift 
clusters including still-forming objects unlikely to be dis- 
covered by x-ray searches. To do this we use the cluster 
red-sequence detection method (Gladdcrs & Yee 2000), 
applied to data from the Red-Sequence Cluster Survey 
(RCS; Gladdcrs & Yee 2005) and the Sloan Digital Sky 
Survey (SDSS, York et al. 2000) Data Release 2 (DR2). 
The RCS data are R c and z' imaging to a depth sufficient 
to detect clusters to z ~ 1.4, more than sufficient to find 
clusters at the redshift considered here. The SDSS data 
are much shallower, but sufficient to find rich/massive 
clusters at the redshift of interest. 

Approximately 50 square degrees of the RCS imaging 
are readily visible from the Las Campanas Observatory 
and we initially searched this area for candidate rich clus- 
ters in the desired 0.3< z <0.5 redshift range. Clusters 
were detected as over-densities on the sky and in color 
and magnitude space. The richest systems were consid- 
ered as candidates for the ICBS program. The lack of 
sufficient RCS imaging area visible from the southern 
hemisphere at 6 h to 9 h and 14 ft to 20 h forced us to use 
the SDSS DR2 data as a secondary source for rich clus- 
ter candidates. In order to ensure a reasonable match in 
mass between the two cluster sub-samples, we restricted 
our attention in the SDSS to areas of low-extinction sky 
comparable in size to the RCS search area, at two RAs 
selected to facilitate the extensive ICBS data collection. 
The total comoving volume covered by this search is 
equal to that in the two Galactic caps out to a redshift 
z ~ 0.055, and should, therefore, include a number of 
rich clusters comparable to those found in local surveys. 
We identified 5 fields which, from the RCS or from our 
analysis of the SDSS data, appeared to contain very rich 
clusters at redshifts between 0.3 and 0.5. Preliminary 
spectroscopy showed that the cluster in one of the five 
fields was not sufficiently rich to be interesting, leaving 
4 fields, two from the RCS and 2 from the SDSS search, 
whose locations are summarized in Table 1. 

2.1. Spectroscopy 
2.1.1. Observations 

Slit masks for each field were populated with objects 
from the photometric catalogs, with a slight preference 
for brighter objects and a strong preference for objects 
brighter than r — 22.5. After the first mask in each field, 
later masks contained unobserved objects plus a num- 
ber of already observed objects, most with poor spectra 
plus some with adequate spectra to use for repeatability 
tests. Most (39 slit masks) spectroscopy was done with 
the full 30' field of the IMACS f/2 camera on the Baade 
Telescope. However, in order to increase the fraction of 
objects observed in the cluster cores, where slit overlap 
issues make it particularly difficult to obtain adequate 
sampling, 3 masks were obtained using the GISMO im- 
age sheer (q.v. Dressier et al. 2011) on IMACS, and 16 
masks were obtained using the LDSS3 spectrograph on 
the Clay Telescope, which has an 8.3' field. Observations 
were done with a mixture of stare mode and nod-and- 
shufilc mode. Typical total integration time per mask 



Table 1 

Intermediate Rcdshift Cluster Fields 



Cluster 


a 


8 


RCS0221 


2 h 21 m 48 s 


-03°46' 


RCS1102 


ll h 02 m 36 s 


-04°40' 


SDSS0845 


8 h 45 m 30 s 


+03° 27' 


SDSS1500 


15 h 00 m 30 s 


+01°53' 



were in the range of 3 to 4 hours, divided into individual 
integrations of 30 to 45 minutes. All IMACS and LDSS3 
spectroscopy, except for the LDP observations discussed 
below, were made at a dispersion of about 2A per pixel, 
resulting in a spectral resolution of about 6A. Typical 
image quality during the spectroscopic and imaging ob- 
servations was 0.6"-0.8". 

The first few slit masks in each field were observed 
with no band limiting filters, providing spectral cover- 
age between 4300A and 9300A. Later observations were 
done with a filter limiting coverage to 4800A to 7800A. 
The shorter spectra allowed more spectra- about 300- 
to be packed onto a single mask, but lost coverage of 
the Ha line at redshifts z > 0.19. In order to obtain 
Ha observations of as many cluster galaxies as possible, 
one mask per field was devoted to observations of already 
discovered cluster members through a 1000 A wide filter 
centered on Ha at the cluster redshift. 

One mask each was obtained in all fields except 
SDSS1500 with the IMACS Low Dispersion Prism 
(LDP). With a mean resolution A/AA ~ 30, the 
LDP spectra cannot be used to measure individual line 
strengths, but are sufficient to obtain redshifts with a 
mean accuracy of about 0.01. Because the spectra are 
very short, of order 1000 objects can be observed on one 
slit mask, to a depth considerably fainter than with grism 
spectroscopy. More information about the LDP prism 
can be found in Kelson et al. (2012). 

The data sets derived from all of these observations 
are referred to in the following as the ICBS data sets. 
For calibration purposes, we also construct a local sam- 
ple of galaxies with optical spectroscopy, and optical 
and 24 fim infrared photometry. This sample consists 
of galaxies with SDSS redshifts between 0.04 and 0.08 
located in the three fields of the SWIRE survey (Lons- 
dale et al. 2003) which are within the SDSS survey area. 
We take 24 /an flux from the SWIRE observations, and 
take spectra, redshifts and optical photometry from the 
SDSS database. 

2.1.2. Spectral Reductions 

With the exception of LDP observations, all IMACS 
and LDSS3 spectra were reduced using the COSMOS 
data package (Oemler et al. 2011). In general, all obser- 
vations of an individual mask, whether from one night 
or from several observing runs were combined into one 
stack of 2-dimensional (flux vs wavelength and slit posi- 
tion) spectral images of individual slits. Using the inter- 
active spectral analysis program viewspectra from the 
COSMOS package, a boxcar 1-dimensional extraction of 
each spectrum was made. Most spectra were extracted 
over a 1" length along the slit, but wider extractions 
were used for particularly diffuse galaxies. LDP spec- 
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tra were reduced using the methods described in Kelson 
et al. (2012). 

Extracted 1-dimensional slit spectra were turned into 
final calibrated spectra in a 3-step process. Firstly, spec- 
tra were put onto a relative flux scale by correcting for the 
instrumental response using observations of spectropho- 
tometric standards. Because the spectral response of 
IMACS is very stable, a mean spectral sensitivity func- 
tion for each instrumental setup was derived by combin- 
ing all standard star spectra taken with this setup. Sec- 
ondly, spectra were corrected for atmospheric absorption 
using the sum of all of the spectra on each slit mask ob- 
tained on a given night. Because the galaxies observed on 
most individual slit masks have a wide range of redshift 
and spectral type, the spectral features of the individual 
galaxies are averaged out in this summed spectra, leaving 
only the spectrum of the atmospheric absorption plus a 
slowly-varying sum of the individual (redshifted) contin- 
uum shapes. This spectrum, after removing the contin- 
uum variations, was divided into each galaxy spectrum to 
remove the atmospheric absorption. Absorption correc- 
tions for spectra observed through one of the Ha filters 
of cluster members could not be performed in this man- 
ner, since the wavelength range was too narrow and be- 
cause all of the objects had similar redshifts. Absorption 
corrections for these masks was performed using sky ab- 
sorption spectra obtained from other observations of the 
same cluster. Thirdly, the spectra are put on an abso- 
lute flux scale by scaling their values using the difference 
between the synthetic r magnitude calculated from each 
spectrum and the total r magnitude of the galaxy. As- 
suming that there are no significant color gradients be- 
tween the roughly 1" square area of the galaxy observed 
with IMACS and the total galaxy, this corrected spec- 
trum represents the total flux of the galaxy. We shall test 
this assumption later. After calibration, spectra of galax- 
ies obtained from multiple slit masks were combined. We 
have found that better results were obtained if spectra 
were combined using a single weight for all data points 
in spectra from each slit mask, rather than using the 
pixel-to-pixel statistical weights obtained as part of the 
reduction process. Relative weights varied by a factor of 
a few between the best and the worst slit masks. The 
final quality of the spectra is, of course, quite variable. 
A small fraction (typically 20%) were too poor to yield a 
redshift. For the remaining, the median signal-to-noise 
ratio, per 2A pixel, at a rest wavelength of 4500A, ranges 
from 30 at r = 19.0 to 7.5 at r = 22.0. 

Redshifts of galaxies were measured using the cross- 
correlation method of Kelson ct al. (2005). Analysis 
of repeat observations of galaxies show a typical error 
<j{z) ~ 2 x 10~ 4 . This error is, to first order, indepen- 
dent of spectral type (absorption or emission lines) or 
galaxy brightness. It is larger than that expected due 
to the typical wavelength errors determined by measur- 
ing atmospheric emission lines cr(z) ~ 6 x 10 -5 , but is 
comparable to that expected due to slit errors. The lo- 
cation of spectra along the slits show a scatter about the 
expected position of about 1 pixel, equivalent to astro- 
metric errors in the galaxy catalogs of 0.2". A scatter of 
1 pixel in the centroid of the galaxies in the dispersion 
direction corresponds to a redshift error <r(z)of3xl0~ 4 . 

The completeness of the RCS1102 redshift sample, typ- 
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Figure 1. Completeness of the RCS1102 sample, black line- 
fraction of catalog with redshifts; red line- fraction of the catalog 
that was observed; green line- fraction of observed objects that 
yielded a redshift. 

ical of all fields, is presented in Figure 1. The total com- 
pleteness, shown by the black curve is the product of 
two factors: the fraction of objects observed, shown by 
the red curve, and the rate of success in obtaining red- 
shifts from the spectra, shown by the green curve. Al- 
though the success rate declines from almost 100% for the 
brightest objects to about 75% by r = 22, most of the 
incompleteness is due to object selection, which varies 
both with magnitude and with position with the field. 
The decline with magnitude in the fraction of objects 
observed, from about 70% for bright objects to 55% at 
r — 22, is completely due to the algorithm which selects 
objects to be included on a slit mask, and which has a 
mild preference for brighter objects. The spatial vari- 
ation in completeness is small. Because of the number 
of masks used in each cluster, and because a concerted 
effort was made to sample well the cluster core, there is 
very little under-sampling of objects in groups and the 
cluster. Typically, what under-sampled regions exist lie 
near the periphery of the field. 

Since it is generally easier to identify redshifts for 
galaxies with strong emission lines, we have checked for 
such a bias in our redshift catalog. To minimize evo- 
lutionary effects (fainter galaxies tend to be at higher 
redshift, and star formation rates increase with redshift- 
q.v. Paper III), we examine the fraction of galaxies with 
EW{[OII\) > 5A in the redshift range 0.15 < z < 0.35. 
Brighter than r = 20.5 absorption spectra increase in fre- 
quency, because early-type galaxies dominate the bright 
end of the luminosity function. However, between r — 
20.5 and r = 23.0, the fraction of emission line galaxies 
is constant, demonstrating that any spectral type bias is 
minimal in our sample. 

Emission and absorption lines in both the ICBS and 
SWIRE samples were measured using viewspectra. 
This is a semi-automatic process, in which gaussian pro- 
files are fit to each line. The wavelength intervals of 
continuum side bands, and of the line itself, are specified 
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in advance. However, one can interact with the fitting 
process to correct for less than perfect fits by altering 
any of the fitting parameters. Output of the fitting in- 
cludes equivalent widths of emission and absorption lines 
as well as fluxes of emission lines. Because the spectra 
have previously been given an absolute calibration, the 
resulting fluxes represent the total line flux for the entire 
galaxy. 

In fitting B aimer emission lines, an attempt was made 
to properly set the continuum in the bottom of the stel- 
lar absorption line; however, given the typical signal-to- 
noise ratios of these spectra, there is a limit to how well 
this can be done. When measuring the H<5 absorption 
line, no attempt was made to remove contamination from 
US emission. Because the H<5 line is often of marginal 
signal-to-noise ratio in these spectra, an improved mea- 
sure of EW(H<5 ) was constructed by combining the 
strength with that of the He line, which lies on top 
of the Ca H line. He was determined from the difference 
between the Ca K line and the Ca H +He equivalent 
width. Empirically it was found that the H<5 equivalent 
width is related to that of H+He and K as: 

H-K < 3 A H5= 1.73 + 0.5(H-K) + O.OA(H-K) 2 
H-K > 3 A H5 = 0.77(H-K) + 1.35 

(1) 

where H—K is the difference, in Angstroms, between 
the equivalent withs oiH+He and K. Our final H<5 values 
are the average of and that derived from H and K. 



T 




EW([OII]) 

Figure 2. Synthetic EW([OII]), as defined by Equations 2 and 
3, versus directly-measured EWQOII]), for objects in the SWIRE 
sample. 

The [Oil] line lies in a rather clean spectral region and 
is easy to measure except for one complication. Some 
galaxies, even some with very strong [Oil], have a very 
weak stellar continuum at 3727A. Slightly incorrect con- 
tinuum levels in the reduced spectra due, for example, to 
sky subtraction errors, can result in quite large fractional 
changes in the continuum value, and therefore quite large 
errors in the equivalent width of [Oil] . To avoid this prob- 
lem, we define a pseudo-continuum at 3700A based on 
an empirical relation between 3700A flux and a combina- 
tion of M(B) magnitude and B-V color, derived from syn- 



thetic photometry of the SWIRE sample of SDSS spec- 
tra, and from this plus measured [Oil] flux, determine 
EW([OII]). The relation between synthetic and directly- 
measured EW([OII]) is presented, for the SWIRE sam- 
ple, in Figure 2; a best fit to this relation is: 

EW{[OII}) = -1.13 x 10- 32 L{OII)/10-° AM37 (2) 

where the 3700A monochromatic magnitude M37 is ap- 
proximated as 

M 37 = Af B +0.20+0.67(5-F-0.09)+0.0359( J B-y-0.09) 2 

(3) 

Analysis of repeat measurements of individual galax- 
ies show that [O II] and H/3 determinations have typical 
errors of 0.07 dex, and US determinations have typical 
errors of 0.13 dex, however these errors are quite depen- 
dent on spectral quality. 
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Figure 3. Synthetic 3" line strengths vs measured 1" line 
strengths, for galaxies in several masks. 

The ICBS spectra are typically extracted from an aper- 
ture 1" square, equivalent to a 5.3 kpc square area in a 
z = 0.40 galaxy. This is considerably smaller than the 
area containing the bulk of the stars of the typical lumi- 
nous galaxy. Since many of the galaxy parameters which 
we derive depend on an extrapolation from the spectro- 
scopically observed area to the total galaxy, any system- 
atic shift between line strengths in the galactic center to 
those of the entire galaxy could lead to systematic errors 
in, for example, star formation rates and internal extinc- 
tion. To test for this, we take 2-dimensional spectra from 



IMACS Cluster Building Survey. I 



5 



several of the best masks and construct synthetic spec- 
tra of an area 3" in diameter. Figure 3 compares 1" and 
3" measures of EW([OII]) and EW(H5 ). The synthetic 
spectra are necessarily quite noisy, resulting in consider- 
able scatter, but it is clear that there in no systematic 
trend of either with area, indicating that our measured 
line strengths are reliable indicators of the total galaxy 
values. 



2.2. Photometry 

Direct imaging in the griz bands was obtained for the 
RCS0221 and SDSS0845 fields with the f/2 camera of 
IMACS. Imaging in the BVRI bands was obtained for 
RCS1102 and SDSS1500 using the Wide Field CCD cam- 
era on the du Pont Telescope. In addition, very deep r- 
band photometry, complete to r = 25.0 was obtained for 
all fields with IMACS. Photometry of the images was per- 
formed using SExTractor (Bertin & Arnauts 1996). Rest 
frame B-V colors and absolute B magnitudes were de- 
rived from the IMACS and du Pont photometry, using k- 
corrections derived from spectral templates constructed 
from the spectra of SDSS galaxies in the SWIRE sample 
described above. Typical color errors at (20.0, 22.0) mag 
are (0.007,0.04) in (r-i), (0.015, 0.07) in (i-z), (0.05, 0.10) 
in (B-V), and (0.04, 0.08) in (V-R). For a small number 
of galaxies with no direct imaging, we have constructed 
synthetic rest frame colors and absolute magnitudes from 
the (fluxed) spectra. 

The RCS0221 and SDSS0845 fields were mapped by 
the Spitzer Space Telescope (Werner et al. 2004) with 
the MIPS instrument (Rieke et al 2004) in the 24 micron 
band. Data were taken in guest observer program 40387, 
PI Dressier. Including overheads, the observations lasted 
11.7 hr per field. For each cluster, the circular 27' di- 
ameter IMACS field of view was tiled with four over- 
lapping MIPS raster-map photometry sequences. This 
covered the area more efficiently than would scan map- 
ping. For almost all the IMAC field of view, at least two 
raster-maps overlap, providing > 980sec exposure time 
per pixel for SDSS0845, and > 1069sec for RCS0221. 

The raw MIPS images were reduced and mosaicked 
using the MIPS Instrument Team Data Analysis Tool 
(DAT) (Gordon et al. 2005). The temporally- varying 
ecliptic foreground was subtracted separately from each 
pointing. Photometry at 24 micron was obtained by fit- 
ting the point spread function (PSF), using the IRAF 
implementation of the DAOPHOT task allstar (Stetson 
1987). The PSF was created empirically from stars in 
each cluster, and an aperture correction was applied as 
in Rigby et al. (2008). 

An extremely bright foreground carbon star at a, S — 
08:45:22, +03:27:09 (J2000) contaminates an area <~ 2.5' 
in radius in the cluster SDSS0845. The SDSS0845 cata- 
logs were edited by hand to remove artifacts caused by 
this star and account for missing survey area obscured 
by the star. Comparing object counts in our fields with 
the deep 24 /xm counts by Papovich et al. (2004) we de- 
termine that our photometry is complete to about 60/xJy, 
equivalent to a star formation rate of about lM j/r _1 at 
z ~ 0.4, with typical errors of about 15/xJy. 



3. DATA ANALYSIS 
3.1. Determination of Galaxy Properties 
3.1.1. Star Formation Rates 

The star formation indicators available to us are the 
optical [OII]3727A, FT3 , and Ha emission lines, and the 
24 (im mid-infrared flux. In principle, the most direct 
measure of star formation rates comes from extinction- 
corrected hydrogen recombination lines (see, e.g. Ken- 
nicutt 1998a for a discussion of the general problem). 
However, our data is inadequate to determine reliable 
extinction corrections. Of those methods available to 
us, many studies have established that using the 24 fim 
flux is the best choice. A number of empirical calibra- 
tions have been made of the correlation of star formation 
rate with infrared luminosity. Since most of the mid-IR 
flux comes from warm dust heated by absorbed UV ra- 
diation from HII regions (see e.g. Wang & Heckman 
1996), one would expect the best correlation to be be- 
tween star formation rate and the IR bolomctric luminos- 
ity. However, that is an observationally difficult quantity, 
and luminosities in either the IRAS 25 band or the 
Spitzer 24 /j,m band are more practical measures. Most 
studies (e.g. Wu et al. 2005, Calzetti et al. 2005, Alonso- 
Hcrrcro et al. 2006, Rieke et al. 2009) have looked at cor- 
relations between bolometric or 24 /an luminosities and 
other measures of star formation rates. Among the most 
sophisticated analyses is that of Rieke et al. (2009), who 
make use of spectral templates to predict 24 (im bolo- 
metric corrections, as a function of IR luminosity and as 
a function of redshift. 

However, as Perez- Gonzalez et al. (2006) have pointed 
out, the true expected correlation should be between IR 
luminosity and the absorbed rather than total UV lumi- 
nosity, since only the former heats the dust. Put an- 
other way, the correlation should be between the star 
formation rate and the sum of the IR luminosity and 
the escaped UV luminosity. Perez-Gonzales et al. (2006), 
Calzetti et al. (2007), Kennicutt et al. (2009, hereafter 
K09) and Calzetti et al. (2010) have all provided calibra- 
tions of this relation. All have substituted the easily- 
observed escaped Ha flux for the unobserved UV flux in 
this analysis. This is not strictly correct, because the ra- 
tio of absorbed to escaped radiation is much higher in the 
UV than it is at Ha . K09 argue that this discrepancy is 
compensated for by other factors. This is not necessar- 
ily true, but the K09 formulation, also used by Calzetti 
et al. (2010) is more convenient than a more strictly cor- 
rect analysis, and the best test of its usefulness is the 
tightness of the correlation between predicted and true 
star formation rates. 

We shall use the relations between 24 fim and 
Ha luminosities and star formation rate given by Equa- 
tion 17 of Calzetti et al. 2010, but with a slight modifica- 
tion to remove the discontinuities in their formulation at 
L(24) = 4 x 10 42 and 1,(24) = 5 x 10 43 . With this mod- 
ification, and renormalizing to a Salpeter IMF, we have 

1,(24) < 4 x 10 42 

SFR = 8.1 x 10~ 42 [L(Ha) + 0.020i(24)] 

4 x 10 42 < L(24) < 5 x 10 43 

SFR = 8.1 x 10- A2 [L(Ha) + 3 x lQ- 9 i(24) 1 - 16 ] 
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L(24) > 5 x 10 43 

SFR = 2.86 x 1CT 43 L(24) (4) 

In this and all following equations for star formation 
rates, luminosities are in units of erg s~ 1 cm~ 2 , star 
formation rates is in units of Mgyr , and we assume a 
mass scale based on a Salpeter IMF. 

At higher redshift the Ha line is often not observable, 
and we have no measurements of it for about two 
thirds of our galaxy sample; [O II] is the best substitute. 
Although the relationship between line strength and 
star formation rate is more straightforward with the 
Balmer lines than with [O II] , the higher-order Balmer 
lines are both weaker than [O II] and also increasingly 
complicated by underlying stellar absorption lines. 
Using [O II] instead of Ha in the Calzetti et al. (2010) 
method should be equally good or better, since the 
extinction at [Oil] is closer to that in the UV. Us- 
ing the data on normal galaxies from Moustakas & 
Kennicutt (2006), as tabulated in Table 2 of K09 
(hereafter called the K09T2 data set), as well as our own 
data, we find best consistency between SFR calculated 
from L(24) and L([OII] ) with that calculated from 
L(24) and L(Ha), with the following set of relations: 

£(24) < 4 x 10 42 

SFR = 8.1 x 10- 42 [1.3£(Oi7) + 0.020L(24)] 

4 x 10 42 < L(24) < 5 x 10 43 

SFR = 8.1 x 10- 42 [1.6L(O77) + 3 x 10~ 9 .L(24) 116 ] 

£(24) > 5 x 10 43 

SFR = 2.86 x 10~ 43 £(24) (5) 



L(24) < 2.5 x 10 42 

log(SFR) = 0.81[log(L24) - 42.40] 

2.5 x 10 42 < £(24) < 5 x 10 43 

log(SFR) = 0.86[/o 5 (£24) - 42.40] (6) 

This relation seems to be valid for galaxies at all redshifts 
and all 24 /im luminosities within our samples. 




1e+41 1e+42 1e+43 



L(24) 

Figure 4. Star formation rate, calculated from Equations 4 and 
5 vs 24 (im luminosity, for galaxies in the K09T2 (red cirlces), 
SWIRE (open circles) and ICBS (black circles) samples. The green 
line is a best fit to the data; the black line represents Equations 4 
and 5 with no optical emission contribution 



The above relations were derived from low-redshift 
galaxies. For higher redshift objects, observed L(24) will 
systematically depart from rest frame L(24), so the appli- 
cation of a k-correction is necessary. Unlike the optical 
case, where one usually has optical colors with which to 
calculate the continuum slopes needed to obtain the k 
correction, mid-IR photometry in adjacent bands is not 
necessarily available. Instead, we use the tabulated k- 
corrections by Rieke et al. (2009) which have been calcu- 
lated from models for star-forming galaxies, which pre- 
dict the continuum shapes and k-corrections as a func- 
tion of star formation rate. 

If no information on optical emission lines is available 
we must fall back on an empirical calibration of SFR 
vs L(24). Figure 4 presents the distribution of the star 
formation rate calculated from Equations 4 or 5 versus 
L(24) for the K09T2, SWIRE and ICBS data sets. The 
black line represents Equations 4 and 5 in the case where 
L(Ha ) and L(OII) and zero, expected when galactic 
extinction is so high that all the UV flux is absorbed and 
reemitted in the IR. As it should, this line follows the 
lower envelope of the galaxy data. The green line rep- 
resents the best fit to the data, and is of the form 



If mid-infrared photometry is not available, the only 
recourse is optical emission lines. In principle, one should 
be able to use line ratios to correct these lines for ex- 
tinction, as has been used to obtain, e.g. L(Ha) to t for 
the K09T2 sample. Unfortunately, this requires better 
data than is usually available for faint galaxies. Exper- 
imenting with the ICBS and even the SWIRE data sets 
demonstrate that using Balmer lines, or Balmer to [Oil] 
ratios to correct the optical line strengths only introduces 
noise, without improving either random or systematic er- 
rors in star formation rates. Instead, we will derive em- 
pirical relations between observed line luminosities and 
star formation rates. 

In Figure 5 we present the relation between star 
formation rate and observed Ha luminosity, for galaxies, 
in the K09T2 sample, and in the SWIRE and ICBS 
samples, where star formation rates were calculated 
using Equations 4 and 5. All three data sets show a 
similar relation between SFR and L24; the best solution 
is shown by the green line. The scatter about the line 
has a dispersion a(log(SSFR)) ~ 0.16 for the K09T2 
sample, and 0.25 larger for the noisier ICBS data. 

L{Ha) < 7 x 10 40 

log(SFR) = l.Q9{log(L(Ha)) - 40.85) 
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L(Ha) > 7 x 10 40 

log(SFR) = lM(log(L(Ha)) - 40.85) 



(7) 



larger for the ICBS data. 

Finally, if even H/?is unobservable, we must fall back 
on [Oil] , for which the calibration is presented in Figure 
7. The best solution for all the data is 




Figure 5. Total star formation rate vs observed Ha luminosity for 
galaxies in the K09T2 sample (red points), the ICBS (black points), 
and SWIRE samples (open circles). The green line is a best fit to 
the data. 



If Ha is not observable, H(3 is the next best choice. In 
Figure 6, we present the relation between star formation 
rate and observed H/3 luminosity, for the same data sets 
as presented in Figure 5; the best solution for the ICBS 
data is. 



log{SFR) = l.l8(log(L(H/3)) - 40.18) 



(8) 
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Figure 6. Total star formation rate vs observed H/3 luminosity for 
galaxies in the K09T2 sample (red points) and in the ICBS (black 
points) and SWIRE samples (open circles). The green line is a best 
fit to the data. 

The scatter in the H/3 determined SFR is larger: 
a{log{SSFR)) ~ 0.25 for the K09T2 sample, somewhat 



log{SFR) = 1.10(log(L{OII)) - 40.39) (9) 

The scatter for [Oil] is a(log{SSFR)) ~ 0.43 for all 
samples. 
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Figure 7. Total star formation rate vs observed [O II] luminosity 
for galaxies in the K09T2 sample (red points) and in the ICBS and 
SWIRE samples (black points). The green line is a best fit to the 
data. 

It must be emphasized that all of these star forma- 
tion rate calibrations, and particularly those using only 
one optical line or IR band, are only claimed to be valid 
for the star formation rates and galaxy types used in 
the calibration. However, the galaxy sets which we use 
should fairly sample normal luminous galaxies are red- 
shifts z < 1.0. The similarity in the relations between, 
for example, the K09T2 and ICBS galaxy samples, which 
contain galaxies at quite different redshifts, observed in 
very different manners, gives us some confidence that 
these calibrations are indeed appropriate for our sam- 
ple, and for any of the other galaxy samples produced by 
surveys of the general galaxy population of the universe. 
They may or may not be equally applicable to unusual 
objects such as ULIRG's, or to extreme dwarf galaxies, 
or to objects at very high redshift. 

3.1.2. Internal Extinction 

For those galaxies with detected 24 fj,m flux, and de- 
tected [Oil] or Ha flux, we can calculate AV em , the ex- 
tinction towards the emission line regions, from the ratio 
the star formation rate calculated from Equations 4 or 5, 
and that obtained from the same equations in the limit 
of £(24) = 0, i.e. the case of zero extinction. How- 
ever, sometimes it will be useful to know the extinction 
towards the total stellar continuum of a galaxy. (For ex- 
ample, in Section 3.2 we shall use the dereddened galaxy 
colors as a starburst criterion). We determine this by 
the method described below, using the SWIRE data set. 
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Figure 8. Observed rest frame B-V colors of galaxies in the 
SWIRE sample vs the quantity M(B)+2.51og(SFR). Large trian- 
gles are objects with A(V) < 0.35. The solid line is the predicted 
relation for delayed exponential models. 

Figure 8 presents the relation, for objects in this sam- 
ple, between the observed rest frame B-V colors and the 
quantity SFRM = M B + 2.5log(SFR), i.e. the star for- 
mation rate per unit blue luminosity in magnitude form. 
Objects with emission line extinctions AV em < 0.35 are 
shown as large triangles. The solid line is the prediction 
of a set of galaxy models computed with the Padova evo- 
lutionary tracks (Bertelli et al. 1994), adopting a Salpcter 
IMF with masses in the range 0.15-120 M Q , and with 
star formation histories following the form introduced by 
Gavazzi et al. (2002): 

SFR ~ 4,e- t2/2r2 (10) 

We shall call these delayed exponential models These 
models use the observed stellar libraries of Jacoby 
ct al. (1984) in the optical (~ 3400 to 7400 A), and they 
were extended to the ultraviolet and infrared with the 
theoretical libraries of Kurucz (1993, private commu- 
nication). They include emission lines formed in HII 
regions, that were calculated using the photoionisation 
code cloudy (Ferland 1996). The nebular component 
was calculated assuming case B recombination, electron 
temperature of 10 4 K and electron density of 10 2 cm~ 3 . 
The source of ionizing photons was assumed to have a 
radius of 15 pc and a mass of 10 4 M Q . 

Galaxies with low extinction lie close to the line; all 
others are redder and fainter (more positive values of 
SFRM) than the line, as would be expected due to ex- 
tinction of the stellar continuum. We assume that each 
galaxy is reddened and dimmed by an amount necessary 
to move it, along a direction parallel to the reddening 
vector, from a location along the line to its current posi- 
tion. We call the V band extinction of the galactic stellar 
population, determined in this way, AV* , and in Figure 9 
present, for the objects in Figure 8, the ratio of A(V) em 
to AV*, vs SFRM. This ratio of the extinction towards 
the stellar population to that towards emission line re- 
gions has been extensively studied by Calzetti, who finds 
an average value of 0.5 (Calzetti 2001). For those galax- 
ies with well determined ratios (i.e. those with significant 
values of A(V), there is a remarkably tight relation be- 



tween AV*/AV em and star formation rate per observed 
luminosity. Galaxies with weak star formation have low 
ratios of AV*/AV em , in other words most of the extinc- 
tion is close to the HII regions, while, in galaxies with 
vigorous star formation the extinction is spread through- 
out the galaxy. 
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Figure 9. ratio of AV, to AV em vs M(B)+2.51og(SFR), for ob- 
jects in Figure 8. open circles- all objects; filled circles- objects with 
AVem > 1.0, which should have more reliable values of AV* /AV em - 
The solid line is the relation summarized in Eq. 11. 

There is a simple explanation of this. Let us make two 
(oversimplified) assumptions: (1) all dust is associated 
with star-forming regions, and (2) the star-forming re- 
gions arc distributed over the volume of the galaxy in 
the same way as the stars. Therefore, the average path 
length for a photon exiting the galaxy will be the same for 
emission line and continuum photons. Now, if / is the 
number of star-forming regions along the line of sight 
of a photon departing the galaxy, then the number of 
dust clouds encountered by a photon from an HII re- 
gion within a star-forming region on the way out of the 
galaxy is P em ~ 1 + /, while the number of dust clouds 
encountered by photons from a star in the general galaxy 
population is ~ /. Thus, AV * /AV em = P*/P em <~ 
//(l + /), which goes from for small / to unity for 
large /. While an undoubted oversimplification, a quali- 
tatively similar trend must exist as long as some fraction 
of the galactic dust is associated with star-forming re- 
gions, which we know is true. 

The solid line in Figure 9 is the relation: 

SFRM < -20 r = 0.12 

-20<SFi2W<-18.1 r = 0.12 + 0.305(SFRM + 20) 
SFRM > -18.1 r = 0.70 

(11) 

3.1.3. Mass determination 

Bell & dc Jong (2001) present simple prescriptions, 
based on Bruzual & Chariot population models, for de- 
termining the mass-to-light ratios of galaxies from op- 
tical broadband colors. As they point out, their pre- 
scription for rest frame B-V colors, for galaxies with 
a Salpeter IMF and solar metallicity: log{M / Lb) — 
—0.51 + 1.45(5 — V), has the great virtue that it is al- 
most parallel to the standard reddening vector. Thus, 
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Table 2 

Polynomial Coefficients for M/Lb vs B-V 



Rcdshift ao ai 02 0,3 0,4 a$ 

0.00-0.25 -4.23 34.64 -121.05 220.29 -197.31 69.03 

0.25-0.55 -1.60 6.85 -9.35 5.22 

0.55-1.00 -1.78 8.14 -12.41 7.28 



although mass-to-light ratios derived by this method 
are very sensitive to galactic extinction, masses are not 
since they are the product of two quantities, L ga i and 
(M/L) ga i with almost exactly inverse dependence of ex- 
tinction. If the predictions of such population models 
are close to correct, our derived masses will be as well. 
The most important dependence in such models is on the 
initial mass function, and we make the same assumption 
in calculating masses- a Salpeter IMF, as we have made 
in calibrating star formation rates. 

We use a variant of the Bell & de Jong approach, but 
calculate mass-to-light ratios using the delayed exponen- 
tial models described in the previous section. We calcu- 
late the predictions at the epochs observed at redshifts 
between 0.0 and 1.1; the results are presented, in com- 
parison with the Bell & de Jong prescription, in Figure 
10. 

Our results are similar to those of Bell & de Jong, but 
they vary with epoch, and a linear relation is not the 
best representation, Instead, we use the polynomials of 

n 

the form log(M/Ls) = ^a,i(B—V) 1 , with values of a, 



presented in Table 2. 

These polynomials do not reflect the full upturn seen 
in the models at the red end. This is deliberate; given 
errors in the observed colors, a relation as steep as the 
model curves would produce much too large values of 
M/L for some galaxies. Since these relations are not 
exactly parallel to the reddening vector, we correct both 
M/L and L ga i for extinction before calculating masses. 
If a continuum extinction values is not available for a 
galaxy, we assume a value A(V)* = 0.4. Because the 
above relations are close to the reddening line, the effect 
of incorrect extinction values on the derived masses will 
be small. 

3.2. Starburst Criteria 

The role that starbursts might play in the evolution 
of galaxy populations, both in clusters and the field, has 
been a subject of considerable attention and controversy, 
at least since Dressier & Gunn (1983). Papers II and III 
of this series (Dressier et al. 2013, Oemler et al. 2013) will 
examine the evidence provided by this survey in some de- 
tail. The ICBS data provide multiple means of detecting 
starbursts. 

We first reexamine the usual spectroscopic indica- 
tors. The equivalent widths of [OII]3727A and US have 
been used for many years as starburst indicators, (e.g. 
Dressier & Gunn 1983, Couch & Sharpies 1987, Dressier 
et al. 1999); too strong an [Oil] line can only be produced 
in galaxies during a starburst, and too strong an H(5 line 
is produced either during or after a starburst. 

We wish to recalibrate the threshold strengths of both 
lines which separate normal from bursting star forma- 
tion, using stellar population models and empirical evi- 




(B-V) 

Figure 10. B band mass-to-light ratios vs (B-V) color. 
The brown line represents the Bell & dejong (2001) solution 
log(M/L B ) = -0.51 + 1.45(B - V). The dashed black lines are 
fits to the delayed exponential model predictions, shown by colored 
points and lines, in each redshift range, a- redshifts of 0.6 (blue), 
0.7 (green), 0.9 (red), and 1.1 (black), the fit is to only the 0.6 and 
0.7 relations, b- redshifts of 0.3 (blue), 0.4 (green), and 0.5 (red), 
c- redshifts of 0.0 (blue), 0.1 (green), and 0.2 (red). 

dence. The behavior of [O II] is fairly simple to under- 
stand. Its equivalent width is the ratio of the emission 
produced by HII regions to the nearby stellar continuum, 
to which stars of all ages contribute, but younger, hot 
stars contribute the most. Dust extinction can only pro- 
duce one possible effect: selectively diminishing the HII 
emission line strength relative to that of the more broadly 
distributed blue stars. Figure 11 presents the EW([OII]) 
vs SSFR distribution of galaxies in the SWIRE and 
K09T2 samples, and of ICBS galaxies. The red lines 
are the predictions of the extinction-free delayed expo- 
nential models for normal galaxies at redshifts (bottom 
to top) of 0.0, 0.5, and 1.0. These predictions can be 
taken as maximum allowed values of EW([OII]), since 
selective extinction of star-forming regions, as described 
in Section 3.1 will generally reduce the observed value of 
EW([OII]) below its extinction-free value. We take as 
the threshold for starbursts the green line, which has the 

n 

form EW([OII}) = £ ai {\og(S S F R))\ with values of a* 


as presented in Table 3. Necessarily an Oil criterion for 
starbursts must miss a significant fraction of them, since 
objects which start off far from the green line may not 
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cross the threshold during even a strong starburst. 
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Table 3 

Polynomial Coefficients for EW([OII]) vs SSFR 



SSFR 


ao 


a i 


02 


«3 


< 3 x 10" 12 


-10 








3 x icr 12 - 4 x icr 10 


3940 


1126 


123.2 


4.033 


> 4 x icr 10 


-56 









value of AEW(H£) , some chance that the object is 
not a starburst, but merely an outlier in the [Oil] - 
US distribution. For dust-free passive galaxies we can 
have more confidence; theoretical models (Poggianti 
et al. 1999) agree with empirical evidence that 3A is an 
upper limit to the US strength of normal galaxies. 
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Figure 11. EW([OII]) versus specific star formation rate. Red 
points: galaxies from SWIRE and K09T2, black points- ICBS 
galaxies. Red lines are the relations predicted by our models for 
redshifts of (top to bottom) 0,0, 0.5, and 1.0. Green line is the 
adapted threshold for starbursts. 

The effect of extinction on the US line is more com- 
plex, and has been studied in some detail by Poggianti 
et al. (2001). In a dust-free system, H<5 will be strongest 
in systems dominated by A stars- such as post-starburst 
galaxies in which the OB stars have died but the A stars 
have not. In systems dominated by older, cooler stars the 
line is weaker, in younger, hotter systems the line begins 
to be filled in by HII region emission. The effect of adding 
dust will depend on the relative distribution of dust, HII 
regions, A stars, and cool stars, and may either enhance 
or diminish the strength of the US line. Some star- 
bursts are known with no visible optical emission lines 
but strong R8 (Smail et al. 1999, Dressier et al. 2009b), 
presumably cases where the HII regions are heavily ab- 
sorbed, but the A star products of the starburst have 
migrated away from the dustiest regions. 

Because of the complexities, a theoretical prediction 
of US strength is impossible for any but dust-free sys- 
tems. For galaxies with significant Population I, an em- 
pirical determination is necessary. US correlates equally 
well with [Oil] strength and with broadband color; we 
choose, following previous practice, to use [O II] . In Fig- 
ure 12 we present the dependence of H<5 on EW([OII]) 
for galaxies in the SWIRE and K09T2 samples. Be- 
cause there is no reliable theoretical prediction, and be- 
cause the distribution is broad, we define a measure 
AEW(H<5) , which is the strength of US relative to the 
line defined in Figure 12, with values, EW([0 II]) <= 
3.0A : EW(RS ) = 3.0; EW([0 II]) > 10A : EW(RS ) = 
5.5, and a linear increase between the two. 

The quantity AEW(H<5) is a measure of the likelihood 
that an object is a starburst; within the SWIRE sample 
about 95% of objects have values of AEW(H<5) less than 
zero. As with our [O II] measure, the AEW(H<5) criterion 
will necessarily miss some fraction of starbursts in 
galaxies whose initial location in the [O II] -H<5 plane 
is far from the AEW(H5) = line. However, un- 
like our Oil measure, there will be, at any positive 
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Figure 12. a- EW(H<5 ) vs EW(JO II]) for galaxies in the SWIRE 
(red points) and ICBS (black points) samples. The green line de- 
fines AEW(H<5) = 0. 

Larson & Tinsley (1978) were the first to demonstrate 
that the color distribution of starburst galaxies is broader 
than that of normal objects. In Figure 13 we plot the dis- 
tribution of rest frame B-V colors, corrected for extinc- 
tion as described in Section 3.1.2 vs SSFR for galaxies in 
the MK06 and SWIRE samples, and, as a green line, our 
delayed exponential model predictions for a redshift of 
0.0. (At higher redshifts, the predictions for high values 
of SSFR move parallel to the z = 0.0 locus, extending the 
line to higher SSFR's and bluer colors.) We have deter- 
mined morphological classifications for as many as possi- 
ble MK06 galaxies, using images from the NASA Extra- 
galactic Database, and have identified AGN's using the 
line strength criteria of Kauffmann et al (2003). Galaxies 
which are morphologically and spectroscopically normal 
are displayed as large black circles, AGN's and morpho- 
logically peculiar galaxies are displayed as red circles. 
SWIRE galaxies, for which no morphology information 
is available, are shown as small blue circles. 

Normal galaxies have a very tight distribution about 
the predicted relation, with a(B—V) corr = 0.035. On the 
other hand, as Larson & Tinsley originally demonstrated, 
peculiar galaxies- almost all of which have morphological 
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peculiarities suggestive of interactions, and therefore oi 
starbursts, scatter more widely. The two black lines, con- 
taining the region of normal galaxies, have the following 
form, where 5 = log{SSFR) + 11: 



5 < : (B — V) 
(B-V) 

5 > : (B — V) 
(B-V) 



,„,„ = 0.69 - 0.1305' 
max = 0.89 - 0.1305 

,„,„ = 0.69-0.3155 
= 0.89-0.3155 



(12) 




Figure 13. The distribution of reddening-corrected rest frame 
B-V colors of galaxies in the MK06 and SWIRE samples, versus 
specific star formation rate, green line- predicted relation. Blue 
points- SWIRE galaxies; black points- normal MK06 galaxies, red 
points- MK06 galaxies with either morphological peculiarities oi 
spectral indicators of AGN activity. The black lines are our chosen 
limits for normal galaxies 

Objects beyond the region defined by these lines are 
presumed to be in some phase of a starburst. Objects 
above the top line, with star formation rates too high for 
their colors, are likely young starbursts. Those below the 
bottom line, with colors too blue for their star formation 
rates, are probably young post-starbursts. Only about 
1 /3 of the morphologically-defined starburst candidates 
lie beyond the normal region. As with the previously 
defined Oil and H<5 starburst indicators, color selection 
can only discover a fraction of starbursts. Since most of 
the outliers have very high specific star formation rates, 
it will be predominantly the stronger starbursts which 
are detected by this criterion. 

3.3. Cluster Properties 

Figure 14 presents redshift-declination pie diagrams for 
galaxies in the 4 fields in the redshift interval 0.10 to 
0.80. The RCS0221 and SDSS0845 fields each contain 
one dominant cluster, at redshifts near where the Red 
Cluster Sequence cluster detection method predicted one 
to lie. The situation in the other two fields is more com- 
plicated; in both there exist two clusters of compara- 
ble richness, and with redshifts similar enough that both 
probably contributed to the RCS detection signal. Pie 
diagrams in the redshift region of each of the 6 major 
clusters are presented in Figure 15. Filled circles repre- 
sent objects that we identify as belonging to the clusters. 
There is, of course, some ambiguity about the correct 
redshift cut in some of these clusters; we have chosen to 
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Figure 14. Redshift-declination pie diagrams for each of the 4 
fields. Objects with redshifts determined by grism specroscopy are 
shown as black points; objects with only LDP spectroscopy are 
shown as open red circles. 



Table 4 

Properties of Clusters 



Cluster 


Ngal 


z 


f20o(Mpc) 


L/L* 


OV200 


Otot 


RCS0221 


245 (337) 


0.431 


1.27 


158 


941 


895 


RCS1102A 


156 


0.255 










RCS1102B 


245 


0.386 


1.39 


224 


772 


697 


SDSS0845 


278 (383) 


0.330 


1.43 


202 


1087 a 


1031 a 


SDSS1500A 


113 


0.420 


1.17 


82 


798 


639 


SDSS1500B 


160 


0.518 


1.23 


182 


844 b 


904 b 



a only includes redshift range 0.317-0.343 
b only includes redshift range 0.507-0.531 

err on the side being too inclusive, so as not to miss any 
potential cluster members in our later analysis, and so as 
not to contaminate our field sample by cluster members. 

With the cluster membership as defined, Table 4 
presents a summary of cluster properties. Numbers in 
parentheses include LDP spectra. The radius r2oo is de- 
fined in the usual way (Carlberg, Yee, & Ellingson 1997), 
a is the velocity dispersion of all cluster members within 
a projected radius of r2oo 7 N ga i is the observed num- 
ber of members, and L/L* is the total cluster luminos- 
ity, in units of L*, calculated assuming a Schechter lu- 
minosity function with the parameters, a = —1.05 and 
Mb* = f(z) as we determine in Paper III, and correct- 
ing for sample incompleteness as a function of magni- 
tude and position. The quantity L/L* should not be 
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over-interpreted: it refers to a volume large compared to 
either the virial radius or the typical radius to which clus- 
ters are normally observed, and therefore includes much 
supercluster population. Values are not given for some 
properties of RCS1102A: the cluster center apparently 
lies close to the edge of the survey area, and we cannot 
know what fraction of the total cluster is observed. 



RCS0221A 




RCS1102A 

:.• * 

. *,y • : 

• m. " $ *: 

°o • V 

• : . 


RCS1102B 

? y . °- 

. • • °ir W 1 ' ' • ° 

•X is l'$ . .• • 
•J. -15 . * • 


SDSS0845A 

° i — ** 

• 




SDSS1500A 

.• • * •:* 


SDSS1500B 

; ft. *. 
*«} *•»«• » v.* 



Figure 15. Rcdshift-dcclination pic diagrams centered on each 
of the 6 massive clusters found in the 4 fields. Objects with only 
LDP redshifts are not included. Open circles are field galaxies; 
filled circles are galaxies assigned to the clusters. The redshift 
range in each plot is z c [ uater ± 0.07, the declination range is 0.7°. 



4. SUMMARY 

We have obtained photometry and usable spectra of 
6002 galaxies in 4 fields of 30' diameter, from which were 



measured absolute magnitudes, rest frame colors, red- 
shifts and absorption and emission line strengths. The 
6002 galaxies includes 1394 members of 5 clusters. Deep 
24 fim Spitzer photometry was also obtained for two of 
the 4 fields. Using new calibrations of star formation 
rates from optical and IR indicators, star formation rates, 
or an upper limit were obtained for 96% of the galaxies 
in the redshift sample (71% detected SFR, 25% upper 
limits). From colors, masses were determined for 87% 
of galaxies with z < 0.7. In addition, at least one mea- 
sure of the presence or absence of an ongoing or recent 
starburst was obtained for 69% of the redshift sample. 

With these data in hand, we shall examine the evo- 
lution of the star formation properties of galaxies in 
the immediately following papers. Paper II (Dressier 
et al. 2013) begins the examination of the processes driv- 
ing the evolution of galaxies infalling into clusters. In 
Paper III (Oemler et al. 2013) we examine star forma- 
tion in field galaxies out to z — 0.6. Paper IV (Glad- 
ders et al. 2013) constructs a more detailed model of field 
galaxy evolution, using the data discussed in Paper III 
as well as the observed evolution of the star formation 
rate density. In forthcoming papers we will elaborate on 
this model, as well as examine the effect of mergers and 
interactions on star formation and evolution, the near 
infrared properties of our galaxy sample, as well as other 
aspects of the cluster environment. 
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